The strength of the radial-breathing mode in single-walled carbon nanotubes 
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We show by ab initio calculations that the electron-phonon coupling matrix element M e - P h of the 
radial breathing mode in single-walled carbon nanotubes depends strongly on tube chirality. For 
nanotubes of the same diameter the coupling strength |-M e - P ft| 2 is up to one order of magnitude 
stronger for zig-zag than for armchair tubes. For (ni,ri2) tubes M e - P h depends on the value of 
(rti — 712) mod 3, which allows to discriminate semiconducting nanotubes with similar diameter 
by their Raman scattering intensity. We show measured resonance Raman profiles of the radial 
breathing mode which support our theoretical predictions. 

PACS numbers: 63.20.Kr, 71.15.Mb, 78.30.Na 



The radial breathing mode (RBM) is without doubt 
the best known feature in the Raman spectra of carbon 
nanotubes. In this vibration all carbon atoms move in 
the radial direction creating a breathing-like deformation 
of the entire tube. This mode is unique to single- walled 
carbon nanotubes and is taken as indicative of the pres- 
ence of nanotubes in a sample. Moreover the frequency of 
the radial breathing mode is proportional to the inverse 
diameter of the tube [jj. Raman scattering is therefore 
often used to determine the diameter or diameter dis- 
tribution in nanotube samples 0, E3- El In detail, the 
relation between nanotube diameters in real samples and 
the radial breathing mode spectrum is more complicated, 
because of the resonances in the Raman process and ad- 
ditional force constants coming from the tube-tube van- 
der-Waals interaction in bundled nanotubes 00,0] ■ Fur- 
thermore, the RBM eigenvector has a small non-radial 
component 0,0. 

It was suggested to use the RBM to find not only the 
tube diameter but also the chiral angle, i.e., to identify 
a particular (774,712) nanotube 0, 0, E3- n i an d 
specify the chiral vector c = ri\a\ +71202 around the cir- 
cumference of a nanotube in terms of the graphene unit 
cell vectors a\ and 0,2. This assignment relied mostly 
on the frequency of the RBM, sometimes combined with 
an argument about the resonant enhancement of the Ra- 
man intensity for the laser excitation energy 0. It was, 
however, always assumed that the electron-phonon cou- 
pling of the RBM is independent of the chirality of a 
tube ■ This means that far from resonance or exactly 
in resonance the scattering intensity of the radial breath- 
ing mode is expected to be the same for different (n\, 77,2) 
nanotubes. Only recently, Strano et al. 0] suspected a 
smaller matrix element for armchair tubes than for zig- 
zag tubes from their measurements of the RBM signal 
strength of a series of carbon nanotubes in solution. 

In this article we show that, contrary to the 
wide-spread assumption, the electron-phonon coupling 
strength of the radial breathing mode depends on the 
diameter and chirality of the nanotube. In ab initio 



calculations we find the squared electron-phonon ma- 
trix elements in zig-zag tubes to be up to one order of 
magnitude stronger than in armchair tubes for the same 
optical transition energy. In semiconducting nanotubes 
the matrix elements allow to distinguish between the 
— 712) mod 3 = ±1 nanotube families. A similar inten- 
sity difference is expected for the two transitions of metal- 
lic nanotubes in each branch of the Kataura plot ^jjj . We 
show experimental evidence of this intensity difference 
based on measurements of resonant Raman profiles of 
the RBM of nanotubes in aqueous solution. The relative 
Raman intensities can independently confirm an (77,1, 772) 
assignment obtained, e.g., by photoluminescence. 

In the expression for the Raman-scattering cross- 
section from perturbation theory the square of the 
electron-phonon matrix elements | A^ e -pft. | 2 appears in the 
numerator. The intensity of the Raman signal in full res- 
onance is scaled by the electron-phonon coupling |l4| . 
When calculating these matrix elements both electrons 
and holes must be taken into account. To every elec- 
tron excited into a conduction band c and interacting 
with a phonon, corresponds a hole in the valence band 
v. Adding up the two contributions, i.e., assuming 
strict electron-hole symmetry, we obtain for the electron- 
phonon matrix element 

Me-ph = (kc\H e - ph \kc) - (kv\H e - ph \kv), (I) 

where the minus sign comes from the opposite charges of 
holes and electrons. 

The diagonal matrix elements of the electron-phonon 
coupling Hamiltonian H e - V h for optical phonons can be 
obtained from the shift of the electronic bands under de- 
formation of the atomic structure corresponding to the 
phonon-pattern 

where the sum runs over all atoms in the unit cell, k and 
b denote, respectively, the wave vector and band index 
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FIG. 1: Calculated band-energy changes (absolute value) per 
unit change in radius for zig-zag (circles) and armchair tubes 
(triangles). Lines are linear fits of the data corresponding 
to the indicated k value (as defined in the horizontal axis of 
Fig-Et); they include the origin as a data point. Inset: Band- 
energy change \dEb(k) /du\ for the first optical transition of 
nanotubes with diameter ~ 8 A as a function of chiral angle. 



of the electronic state, i indexes the phonon, M is the 
atomic mass, e l a the polarization vector of the phonon, 
normalized as ^ a e l a e J a = r5y , -Efc(k) the electronic energy 
and u a the atomic displacement. N is the number of unit 
cells in the system (1 in our calculation). 

We calculated the T-point phonon spectrum and the 
band structure of several isolated nanotubes in their 
minimum-energy configuration and under deformation 
due to the RBM to obtain the change of the electronic 
energies in Eq. © 16]. The RBM eigenvector was ob- 
tained from finite differences; it has a small non-radial 
component 0, II] which lowers the calculated matrix ele- 
ments by « 30 %. All calculations were performed with 
the SIESTA code 0] within the local density approxi- 
mation |l8| . The core electrons were replaced by non- 
local norm-conserving pseudopotentials [lj|. A double- 
C, singly polarized basis set of localized atomic orbitals 
was used for the valence electrons, with cutoff radii of 
5.12 a. u. for the s and 6.25 a. u. for the p and d or- 
bitals 20]. 16 k points in the k z direction were included 
for metallic nanotubes and 3 k points for semiconducting 
tubes. Real-space integrations were done in a grid with 
a cutoff w 270 Ry. 

In Fig. |I]we show the energy change dEt,(k)/du = 
d[E£(k) — E£(k)]/du for several tubes as a function of 
inverse tube diameter. The data was evaluated at the re- 
gion of the Brillouin zone with the highest optical transi- 
tion probability, i.e., the band extrema HU • We grouped 
the data corresponding to k points that are close when 
mapped onto the graphenc Brillouin zone (lines in Fig.^l. 
To first approximation, the points of a particular group 
correspond to the same transition energies. We find that 
for a particular transition energy, \dEb(k)/du\ is propor- 
tional to 1/d and tends to zero in the infinite-diameter 
limit. This trend can be easily understood, since the 
same change in radius yields smaller bond- length changes 



for bigger tube radii. The infinite-diameter limit corre- 
sponds to a translation of graphene, which cannot affect 
the electronic system. It is clearly seen that \dEt,(k)/du\ 
for zig-zag tubes is up to 2.5 times larger than for arm- 
chair tubes. The inset of Fig. ^ shows the energy change 
for the first transition of tubes all with diameter «8 A 
but different chiralities; in particular, the (8,4) tube lies 
between the armchair and the zig-zag values. The as- 
sumption of a chirality-independent electron-phonon in- 
teraction is thus incorrect. Our results suggest the use of 
relative Raman intensities for discriminating chiralities. 

In Table HI the calculated matrix elements and RBM 
frequencies are summarized; we found wrbm ~ Ci/d+C^ 
with C\ = 232 cm _1 nm in excellent agreement with the 
literature (H and C 2 = -6 cm' 1 . The largest dif- 
ference in the matrix elements A4 e - P h between zig-zag 
(bold face) and armchair nanotubes of similar diameter 
and, hence, same RBM frequency is found for the (11,11) 
and the (19,0) tubes. The matrix element .M3 of the 
(19,0) tube is by a factor of three larger than .Mi of 
the (11,11) tube although the two transition energies are 
similar (see Table [IJ. Since the Raman signal is propor- 
tional to |A^ e -p/i| 2 , we expect the RBM intensity to be 
nine times larger for the (19,0) than for the (11,11) nan- 
otube from the difference in |-M e - p /j| alone. For different 
transition energies this ratio could be even larger. 

The matrix elements of zig-zag tubes show another in- 
teresting feature: they have either a larger magnitude 
and are negative or a smaller magnitude and are positive 
[e.g., for the (10,0) tube Mi = -0.028 and M 2 =+G.017]. 
A change in sign is very uncommon in electron-phonon in- 
teraction in solid-state systems. The matrix elements are 
positive in most semiconductors [24| . To explain this un- 
usual behavior we calculated dE^ 1 '(k) / du of a graphenc 
sheet (aQ qulL = 2.467 A) stretching the sheet in the zig- 
zag direction to simulate the radial atomic displacement 
and adding the non-radial component (see Fig EJ) . 

In Fig. |2t we show \dEb(\t) / du\ for graphene when 
stretching it according to a (19,0) tube (solid line) and 
a (17,0) tube (dotted) together with the ab initio calcu- 
lated values for these two tubes. M. e -ph depends in sign 
and magnitude on where the optical transition occurs 
with respect to the if -point of graphene. The r-point 
states of an (n, 0) zig-zag tube in the graphene BZ are 
obtained by dividing the T-K-M line into n parts (see 
Fig. ■ The states closest to the if -point of graphene 
have the lowest transition energies. The (17,0) tube, e.g., 
has its 1 st transition to the left of the if -point and the 
2 nd one to the right. The energy shift of the graphene 
bands is smaller to the left of the K point than to its 
right (Fig. Et)- Therefore for this tube \Mi/M 2 \ < 1. 
The (19,0) tube, on the other hand, has its I s * transition 
to the right of the if -point and the 2 nd to the left, yield- 
ing \M.i/Ai2\ > 1- Furthermore, dEt(k)/du is negative 
to the right of the if -point, and positive to the left, ex- 
plaining the signs of A4 e -ph- In general, all semiconduct- 
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446 


287 
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179 


170 


151 
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Mi 


0.050 


-0.028 


-0.015 


-0.013 


0.021 


-0.010 


0.016 


-0.022 


-0.017 


0.014 


-0.005 


-0.015 
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(0.7) 
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(0.6) 


(1.4) 


(0.5) 


M 2 
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0.017 
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(2.5) 
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M 3 
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-0.017 
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-0.017 


-0.012 
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(2.6) 


(2.6) 
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(2.6) 


(1.9) 


(2.1) 


(3.3) 


(1.6) 


Ma 




-0.031 

(3.0) 






-0.028 
(3.1) 






-0.022 

(3.2) 




0.009 

(2.4) 




-0.016 

(2.7) 



TABLE I: Calculated diameters, RBM frequencies, and electron-phonon matrix elements M e - P h (in eV) for the first optical 
transitions (ab initio calculated energies in eV in parenthesis The matrix elements (bold for zig-zag tubes) were rounded 

to 0.001 eV. Rows labeled Mi-i correspond to the first four optical transitions for light polarized parallel to the nanotube axis. 
A '-' indicates a lack of the transition or a band shift which could not be evaluated for technical reasons. Note that M e - P h is 
normalized to the unit-cell volume V c . 
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FIG. 2: a) Detail of the Brillouin zone (BZ) of graphene 
around K with the band lines of a (19,0) nanotube and the 
Y-K-M line, b) Electronic bands of graphene along T-K-M 
in equilibrium (solid line) and under a deformation of 0.1 A 
corresponding to the RBM of a (19,0) tube (dashed, differ- 
ence enhanced xlO). c) Calculated dEt(k) /du for a (17,0) 
(squares) and a (19,0) (circles) tube, and of graphene de- 
formed to simulate the RBM of a (19,0) tube (solid line) and 
a (17,0) nanotube (dotted). 



ing tubes can be divided into v = (ni — n 2 ) mod 3 = ±1 
families, which behave like the (17,0) and (19,0) tube 
with respect to sign and relative magnitude of the A4 e - P h- 
Metallic nanotubes usually have two close-by transition 
energies due to trigonal warping |25| . The transition with 
lower energy originates from the right of the if -point, 
the one with higher energy from its left. Therefore, the 
lower-energy transition is expected to give a higher Ra- 
man intensity. 

To confirm our theoretical predictions we performed 
Raman scattering measurements on nanotubes in solu- 
tion [2rl |27| . Raman spectra were excited with a Ti- 
Sapphire laser, recorded with a DILOR XY800 spectrom- 



eter, and corrected for the sensitivity of the experimental 
setup. We then calculated the squared scattering ampli- 
tudes |W^/| 2 from the Raman signal. 

In Fig. [3^ we show two selected resonance profiles of 
radial breathing modes. Using the assignment by Bachilo 
et al. [2(| we identify these resonances as the second tran- 
sition of the (14,1) nanotube with v = +1 and the (11,0) 
tube with v = — 1 |2j|. We predicted a higher Raman 
intensity for nanotubes with v = —1 in excellent agree- 
ment with the experimental data. As shown in Fig. ^ 
dEb(k)/du is similar for nanotubes with similar diame- 
ter, chiral angle, and v. Approximating the matrix el- 
ement for the (14,1) tube by dE b (k)/du for the (16,0) 

i 3. 



nanotube we find theoretically | / 



Experimentally, the ratio {Wpj'^/Wpj' 1 ^ \ 2 ~ 4 is in ex- 
cellent agreement with our ab-initio result and a uniform 
distribution of chiral angles in nanotube samples [2j| . We 
stress that the intensity difference between v = — 1 and 
+1 tubes is generally observed in our experiment and not 
limited to the two profiles shown. A more detailed study 
is underway, but beyond the scope of this paper. 

Raman-based (ni, n 2 ) assignments performed so far re- 
lied strongly on a Raman intensity analysis H that im- 
plicitly assumed constant electron-phonon coupling. As 
shown, this is not correct. The most intense peak does 
not necessarily correspond to the nanotube closest to res- 
onance. As an example, we show in Fig. |2b a Raman 
spectrum for an excitation energy of 1.65 eV. This exci- 
tation energy corresponds to the maximum in the reso- 
nance profile of the (14,1) nanotube, whereas the (11,0) 
nanotube is excited 25 meV below its resonance. The 
Raman intensity thus depends not only on the resonance 
condition, but also on the particular tubes under study. 
Differences in Raman intensity due to resonances can- 
not be distinguished from the chiral angle dependence of 
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FIG. 3: a) Measured resonance Raman profiles (symbols) and 
fits (lines) of two different nanotubes. b) Raman spectrum 
with a laser energy of 1.65 eV. 



M-e-ph using only a single excitation energy; instead a 
resonance profile has to be evaluated. 

A chirality-dependent electron-phonon coupling natu- 
rally explains the observations by Strano et al. 12] . They 
mapped the electronic transitions of metallic nanotubes 
using Raman excitation profiles. Strangely, the armchair 
tubes were apparently missing in their sample. In con- 
trast, photoluminescence experiments on the same type 
of sample [26| suggested a predominance of large-chiral- 
angle tubes. Our calculations solve this apparent con- 
tradiction: the Raman signal of armchair nanotubes is 
small due to a weak electron-phonon coupling. 

In conclusion, we calculated the electron-phonon ma- 
trix elements for carbon nanotubes. The matrix ele- 
ments of zig-zag tubes are much larger than those of 
armchair tubes, leading to a larger Raman signal for 
smaller chiral angle tubes. Furthermore, for semicon- 
ducting tubes the magnitude and the sign of the ma- 
trix elements change systematically for different optical 
transitions and v = ±1 nanotubes. Relative Raman in- 
tensities can discriminate between armchair and zig-zag 
tubes as well as v = ±1 tube families. The latter we 
demonstrated by measuring the radial breathing mode 
resonance on a —1 and +1 nanotube. The family and 
chiral-angle dependence of the Raman intensities can be 
used for a refined assignment of chiral indices and chiral- 
ity distributions. 
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